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Abstract 

We study the complex-valued resonance spectrum of a dc-SQUID coupled to a flux qubit, where 
the former is treated in the cubic and the latter in the two-level approximation. It is shown 
that this spectrum is well-defined and contains most of the relevant information on the escape 
process. Thus, the language of resonance states is precise and well-adapted to switching- (or 
trigger-) type qubit readout, and a worthwhile complement to the various descriptions of continuous 
qubit measurement. Initial progress is analytic, but nonperturbative numerical methods have been 
formulated and should soon yield accurate results for all parameter values. 
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I. INTRODUCTION 



Parallel to, and stimulated by, the remarkable development of experimental quantum 
computing, there has been considerable theoretical attention for the problem of measuring 
a single two-level system (qubit, pseudospin). There is considerable variety both in the 
studied detector candidates (quantum point contacts [I, Q, 0, 0| , normal 1^ and super- 
conducting single-electron transistors in various regimes) and in the employed theoretical 
apparatus (correlation functions and diagrammatic methods [4, 6], generalized Master equa- 
tions [H, 1^, quantum trajectories j^, and Bayesian formalism [2j). However, the theoretical 
emphasis seems to be on continuous readout while experimentally, trigger-type detectors are 
dominant, at least until recently. In the latter, exemplified by switching-current measure- 
ment in a dc-SQUID, the readout is a single event instead of a continuous process, making 
the theoretical description rather different. Yet, many physical questions, such as how to 
characterize detector efficiency, and what exactly is measured when reading out a super- 
position state in the presence of a non-negligible qubit Hamiltonian, are similar. It is our 
purpose to study these issues, using a language especially suited to trigger- type readout. 

The model consists of a dc-SQUID coupled to a flux qubit, described by the Hamiltonian 

H = --^^ + a<p- P<f - i(ea,+Aa,) + \a,<P . (1) 

The first three terms represent a small-inductance dc-SQUID for a near-critical external 
current bias, where the metastable wells of the tilted-washboard potential are widely spaced 
relative to their width and depth, so that for the escape process one can focus on a single 
such well. Clearly, one needs C > for the capacitance in order to have a well-behaved 
mass term; taking /3 > establishes — > oo as the unstable direction, upon which one 
also needs a > to have a potential well near = 0. We also take A > for the qubit 
tunneling amplitude, and A > for the SQUID-qubit coupling. The latter breaks the 
equivalence of cx^ = ±1, so it is better to let e have either sign. As written, the model (jlj 
has six parameters,^ but it will be shown in Sec. Ill Al that there are only four nontrivial 
dimensionless ones. 

A remark is in place concerning the choice of coupling term Acr^^: while represents 
the qubit flux in the two-state approximation, in (P) is the overall phase difference in 
the dc-SQUID (conventionally, the average of the two junction phases), not the latter's 
loop current (related to the difference of these phases), which in fact gets eliminated when 
modeling a small-inductance SQUID with a single degree of freedom. Thus, the term does 
not simply represent inductive fiux-fiux coupling. One can, however, reason as follows: while 
(3 is fixed for a given SQUID, a is proportional to the (small) difference between the bias 
and effective critical currents, where the latter is dependent on the external flux applied to 
the SQUID; in its turn, this flux may effectively contain a component inductively coupled 
into the SQUID loop by the qubit circuit. If the external influence on a is linearized, this 
mechanism is reflected by our choice (a + \az)(t> in (HJ. Foremost, though, the model (HJ 
is meant to be physically motivated and well-defined for all parameter values; note that a 
metastable and essentially irreversible trigger, coupled to a quantum two-level system, may 
also be a prototype for, say, particle detectors. For an accurate description of a particular 



Achieved in part by taking the inflection point of the potential as the origin of (j), which eliminates the 0' 
and (jy^ terms. 
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(strongly coupled) SQUID-qubit system, it is advisable to do a more systematic derivation, 
starting from a multi-loop, multi-junction Hamiltonian. 

The SQUID potential U{(j)) = acp — allows for outgoing states which never return to 
the well, and a well-defined escape problem results without need for a heat bath or other 
model of dissipation.^ In this setting, one expects the existence of exact resonance states or 
quasinormal modes 

i*(^-*)>=(^:;(^:t))=i*'*"'"'"'' 

where uj is complex with \m.uj < 0. Thus, the exponential in Q accounts for both an 
"energy" Reu; and a decay rate |Ima;| (for the amplitude; the probability decay rate is 
twice this). The notation is merely shorthand for the two-component wave function; 
in particular, there is no implication that is normalizable, although some linear-space 
structure will be introduced in Sec. IIIBI This also means that, while and u should 
follow from the two-component eigenvalue equation 

H\-^) = uj\^), (3) 

the latter's familiarity is in part deceptive. While we thus deviate from textbook Schrodinger 
theory, exact resonances are known for, e.g., Fokker-Planck (generalized diffusion) prob- 
lems jsl, in cavity optics, and in linearized black- hole gravity j9|. If one has an effective 
way to find these resonances also for Q, the above defines a framework for switching-type 
(discontinuous) quantum measurement of perhaps unprecedented rigour and precision, for 
arbitrarily large escape rates. An elementary introduction to quasinormal modes of the 
Schrodinger equation is given in the appendix. 

The challenge becomes to calculate the resonance spectrum {ujj{C,a, P,e, A, X)}, in a 
setting where the standard Schrodinger toolbox is called into doubt. Based on the large-0 
asymptotics, it is suggested to study © not along the real 0-axis, but along a contour in 
the complex plane for which the boundary values are more readily implemented. In Sec. IHl 
this is taken up for the noninteracting 

Hs^-'4^^a^-P^r (4) 

This contour approach is related to, and hopefully a refinement of, the "complex scaling" 
method in ^3|; at the least, results will be compared. Since the problem's instability/non- 
hermiticity is solely due to the SQUID, (jH) should capture the essence of these aspects, upon 
which addition of the qubit presents only technical complications, addressed in Sec. IIII CI 
(Of course, the qubit complicates the interpretation of the results.) 

With the uncoupled problem under control. Sec. IIIII returns to the interacting Hamil- 
tonian (HI). In Sees. IIII Al and IIIIBl this is first of all done perturbatively; interestingly, 
the analysis applies a bit beyond the weak-coupling regime, in that a finite quantum non- 
demolition (QND) qubit-detector interaction is admissible. Subsequently, in Sec. IIII CI a 



^ Accurately modeling a given experimental situation may still require the inclusion of such a heat bath, 
of course. In particular, a switching-current distribution wider than predicted by may indicate the 
presence of additional noise sources. 
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numerical scheme for the interacting case is formulated, the results of which should pro- 
vide general and accurate information, and allow comparison with analytical results such as 



m 



As alluded to above, this work is motivated by questions about the interpretation of 
switching-current measurement. The qubit will affect the escape process in the SQUID so 
the latter may be said to read out the former, but how exactly? A plausible scenario, for 
which there is some experimental support is the "self-averaging" one: for each qubit 
eigenstate k, the SQUID sees the expectation of the qubit flux in that particular eigenstate, 
i.e., a (— > a + {az)f.X. The underlying physical picture is that the influence of the qubit on the 
SQUID is to be averaged over many coherent oscillations in the former. However, one can also 
argue in favour of a "weakest- link" scenario , in which the SQUID predominantly escapes 
when the qubit is in the state o"^ = —1, facilitating said escape, i.e., a i— > a — A whenever the 
state (Tj; = — 1 has an appreciable probability of being occupied. While less prominent in the 
literature, this is the behaviour one would expect of a readout SQUID coupled to a classical 
oscillator. At least for a linear detector, in Sec. IIII Bl these two scenarios are recognized 
as opposite limits of a more comprehensive description. Besides this practical relevance for 
SQUID-based flux measurement, the aim is to contribute to quantum measurement theory 
in general, and show how the same concepts surface in radically different readout methods. 



II. FREE SQUID 

A. Preliminaries and complex-contour formulation 

Consider {Hs — uj)il){(l)) = 0, with ifg as in (@)). The full titled-washboard Hamiltonian 
has three parameters. Eliminating the overall energy scale, one would be left with two 
dimensionless parameters: one characterizing how close one is to critical bias (a purely 
classical property), and one counting the number of levels in each well (also involving the 
capacitance). In the cubic approximation, one always is infinitesimally close to critical bias, 
so one expects only the latter parameter to remain.^ Indeed, scaling s = K,(f), one finds that 
for K = {P'^y^^ (with 7 = C/2e^) the Schrodinger equation takes the form 

-'ip" + {as -s^-u)'ip = , (5) 

where now of course ip" = and where the new parameters are related to the old ones as 
a = (7^/^//3^/^)a and u = {'j^^^ / P'^^^)u. While one could also have left the sole parameter 
in front of either the kinetic or the cubic term, the form has two slight advantages: (A) it 
still has the flexibility to tune the potential from metastable via critical to unstable, and 
(B) the leading asymptotic behaviour for |s| — > oo is independent of a, see around © below. 
Since there is no need to refer to the unsealed problem any more except for comparing to 
experimental parameters, we will proceed using the form (0) and henceforth omit the tildes. 

Next, let us study (0) asymptotically. The leading WKB behaviour is ip^s) ~ 
exp[± J'^du ^JU{u) — uj\ . For each cj, define ipi^^s^uj) and '?/'r(s, cu), satisfying the "left" 
(s — i> — oo) and "right" (s — >• oo) boundary conditions, respectively, but in general not both. 
For t/'l, the potential U{s) oo presents a conventional wall, and one clearly has to take 



^ One can also consider the cubic problem as a harmonic oscillator with a cubic "correction" . Since all 
harmonic oscillators are created equal, only the magnitude of the cubic term is a nontrivial parameter. 
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the bounded ^'his) ~ e 5( «)^^^_ Yot the "right" solution Res — oo, matters are more 
comphcated, since neither solution 



^r(s) ~ exp 



±i / du ^/uj — U{u) 



(6) 



dominates the other [i.e., the positive real s-axis is an anti-Stokes line of (jS))]. One now 
reasons as follows: since ipni^y^) G{s,u;uj) for s > m [see (IT^ below] and G is the 
retarded Green function, 7/^^(5,0;) should be bounded in the upper half tu-plane. Hence, 
one should take the plus sign^ in (jUI), i.e., iPk{s) ~ e*^*^''^. Further, one now finds that 
'0r(s) ~ e^sl^l'' ^ along the line args = vr/S, the Stokes line in the upper-half s-plane closest 
to s > 0. Hence, it is suggested that any numerical analysis proceed asymptotically along 
this line, for maximum stability and minimum oscillation of the solution. Because of the 
general analyticity properties of solutions to linear analytic ODEs, and the nearly ideally 
simple behaviour of U{s) in the complex s-plane, the choice of contour C joining the "left" 
and "right" asymptotes is now flexible and can be based on numerical convenience. A 
hyperbola comes to mind, but two straight lines joining at the origin may be even simpler. 

The next step is actually solving Hsip = ^i'- It is stressed that the above WKB asymp- 
totics were intended only for finding the contour C; the eigenvalues themselves are to be 
found numerically, with in principle arbitrary precision. At least two methods seem viable. 
(A) Discretization: the continuous problem (0) is converted into a matrix problem, which 
is subsequently diagonalized either fully or partially. Care has to be taken when covering 
the infinite s-contour with a finite number of grid points. (B) Searching for roots of the 
Wronskian 

which has the standard property of being independent of the choice of matching point Sq.^ 
At the zeroes of W{uj), which can be found with a standard complex-root solver,^ ipi^ and V'r 
are linearly dependent, which means that either function satisfies both boundary conditions, 
i.e., u is an eigenvalue. An advantage of (B) is that the numerical determinations of V'l and 
ijj-R may be optimized separately. The Wronskian method is also particularly suitable when 
only a few resonances are needed, and their approximate location is known. This is the case 
here, for one is only interested in the few "physical" resonances corresponding to the finite 
number of metastable well states.^ 

Let us finally point out that formulating the problem on a Stokes contour not only makes 
'^L,R well-defined, but also renders their numerical calculation very accurate and stable. 



^ This selection of the outgoing wave also follows by requiring that, for real oj and s — s- +00 in one 
should have J > for the probability current J = — ip*Vip). 

^ For an s-contour consisting of two straight lines, the point where they join is the obvious choice. 

^ Minimization of should be a feasible elementary strategy: not only does \W\'^ clearly attain a 

minimum at zeroes of W, but the standard maximum/minimum modulus theorem of complex analysis 
(readily verified using a Taylor expansion of W) guarantees that these are the only minima. 

^ There is, however, a theoretical interest in the fate of these resonances as they "fall out of the well" when 
a is lowered. Note that the resonance frequencies, being complex zeroes of the analytical function W{uj), 
should in all cases depend continuously on the system parameters according to the argument principle. 
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Namely, choosing initial points si with ^ 1, any error in the initial conditions® amounts 
to an admixture of the undesired dominant solution. The latter is exponentially suppressed 
when integrating "inward" from the asymptotic point si to the matching point sq, while the 
desired subdominant solution is amplified by a similar factor. The appendix illustrates this 
in an elementary setting. 



B. Analytics 

For the potential U{s) = as — s^, the well minimum is at = —^/a/3. From f/"(sw) = 
^cUp one finds Up = 2(3a)^/^ for the plasma frequency. By inversion symmetry, the barrier 
height is AU = — 2f/(s„) = |v3q^. For the number of states in the well, one now has the 

estimate = AU/cUp = |v3a^; if A^s ^ 1? one is within the semiclassical regime, at least 
for the lowest few states n = 0, 1, 2, . . .. In this case, one can refer to the literature for a 
prediction for the resonance eigenvalues: 

oon ~ U{s^) + (n+i)a;p + oof' - f r„(^p, AU) . (8) 

Here, o;^"^^ is a small (in the semiclassical limit) correction due to the anharmonicity of the 
well, the lowest few of which are known explicitly: uJq^^ = — lldj"^, tuj^"^^ = — Tlu;"^, and 
^anh ^ _i9i^-4 xhe decay rate is [3,E3 



(9) 



27m! 

it occurs in (jSI) with a prefactor ^ due to the conversion from probabilities to amplitudes. 
Combining the above, the prediction (jS)) is seen to be completely explicit in terms of a; it 
should be valid whenever n Ng. 

The asymptotics around (jH)) have an interesting consequence: iphis) ~ e^t(~'*)^''^ can be 
extended to |7r < args < |7r, and iI>r{s) ~ e^t^'^^^ to — |7r < args < |7r. Hence, there is a 
region of overlap, and for an eigenfunction ipnis), we conclude that the "left" and "right" 
asymptotic forms should occur with the same prefactor. Note that, for |7r < args < |7r, 
we have compared two dominant asymptotic solutions, which can be done only up to a 
subdominant solution. Therefore, we first have to ascertain by other means that a solution 
satisfying both boundary conditions actually exists, and asymptotics alone cannot yield the 
eigenvalues cUn- In particular, if a solution, say ^pl^{s^—oo,uJ) ~ 1 ■ e~i*^~**^^''^, is extended 

■25/2 '25/2 

to s > for G M, i.e., ?/'l(s^+cxd, cu) ~ ae*5* + a*e~^5'^ , then it seems that the 
preceding argument cannot determine a, because the last asymptotic expansion is valid only 
on |args| < |7r due to the Stokes phenomenon. Of course, these claims should be verified 
numerically. 

For an alternative approach to (jH)), consider the Riccati formulation. If ip" = [U — u}]iJj, 
then 

f' + f = U-u (10) 



From © and above, V'l(*i) = (— si)^/^i/)l(si) and ipji{si) = isi^'^ipB.{si) are seen to be suitable starting 
values. 
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for the logarithmic derivative f = if)'/ if). Thus, a second-order hnear equation has been 
converted to a first-order nonhnear one. This has several advantages. First, the hnearity 
of (0), while a blessing in many contexts, also means that solutions are arbitrary up to an 
overall constant; this ambiguity cancels in the definition of /. For instance, the Riccati 
version of the matching criterion reads simply /l = /r. Also, the WKB approximation 
is obvious by neglecting /' in (llUj) . More importantly, the exponentially large or small 
tails typically present in ip cancel in /, so that one deals with at most algebraically large 
quantities. A severe disadvantage of (fTUj) usually is that it converts the nodes of ip into 
simple poles of /. However, while nodes are ubiquitous for the real ip one usually encounters 
for real s and u;, presently these nodes are readily circumvented in the s-plane. 

The above leads to slightly different recommendations for the integration contour C in 
the Schrodinger and Riccati formulations. For (jSJ, one could either leave the real s-axis at 
s = (upon which yet another reparametrization s = —re^'^^^ reduces the inward problem 
for ^/'r to standard form) to keep the parametrization as simple as possible, or stay on this 
axis at least until the barrier point s = |s„|, so that the (complex) resonance wave function 
is as amenable as possible to a conventional interpretation. For (jlUj) . on the other hand, one 
should leave the real s-axis at an s < to the left of the classically allowed region, where 
will have (near-)nodes on the real axis. 

Some familiar properties of quantum mechanics can be generalized to the present setting. 
Consider wave functions ip and X; each satisfying both boundary conditions on C, and in- 
troduce their bilinear product (V', x) = Jc d'S4'{s)x{s). One now derives (V^, Hx) = (x? Hip), 
through the usual integration by parts. In particular, for eigenfunctions this implies 

{uJm-UJn){i'm,1pn) = 0, Or 

{ipm,ipn) = 0, m^n, (11) 

where we also used that the ID Schrodinger equation has no degeneracies, so m 7^ n if 
uJm 7^ i^n- The derivation of this "orthogonality" relied on the symmetry of H; for the 
latter, it is essential that the bilinear product does not involve complex conjugation, since 
H is not real on C. However, this has the consequence that the diagonal entries {ip,ip) are 
not positive definite, not even necessarily real or nonzero which is why we refrained 
from normalizing "(V'm,'^n) = ^mn" ■ Also, there is no suggestion of completeness. 

For future reference we introduce the Green function G{s, u; uj) = G{u, s; u) through 

N ipL{s,oj)^R{u,u;)e{u-s) + ipi^{u,uj)ipR{s,uj)e{s-u) 

G{s, u; UJ) = — — [s, M G M) , (12) 

W{uj) 

solving [-f^s('5) —uj]G{s, u] lj) = —6{s—u). Since W{u) has zeroes at the resonances a;„, these 
correspond to poles of G, and the eigenfunctions ipn are the associated residues A major 
feature of this work is the generalization to complex s, so that the definition of G potentially 
leads to functions of two complex variables. This is emphatically avoided, and below we 
will at most consider G{s, u; uj) with s and u on the same fixed contour C (having a real 
parametrization), in which case the generalization of (jl2j) is still obvious. 

For an application of the above, let us consider the response to changes in a. This provides 
an opportunity to acquaint oneself with the perturbation theory required in the interacting 
case, but in a scalar (as opposed to spinor) setting. Thus, consider the a-derivative of (0), 

-C + {as -s^- uj)^ + (s - uj')ip = , (13) 

where ^(s, a) = daij{s, a) and uj'{a) = duj{a)/da without risk of confusion, and where ^, ip, 
and UJ still have an implicit dependence on the discrete mode index n. Clearly, defines ^ 
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at most up to a multiple of ifj, reflecting the freedom to choose an a-dependent normalization 
constant in ■?/'(«). Rewriting (fT^ as [{^/ipYip^]' = {s — uj')iIj'^, one sees that the integral of 
the Ihs must vanish over the integration contour C, since ip obeys nodal conditions at both 
ends for any a, and therefore so does ^. Hence, the same must hold for the integral of the 
rhs, i.e., 

duj _ J^dssij^s) _ {ij,stlj) _ 

da ysij^is) " ■ ^ ' 

This relation gives the a-response in a familiar form, as the normalized matrix element of 
the perturbing operator s. However, the integrals are evaluated on the contour C joining 
two Stokes lines, ensuring rapid convergence, and importantly involve tp"^ not {ipl"^- Note also 
that ()14|) simultaneously gives the real and imaginary parts of the shift. 



C. Numerical results 

So far, the integration of '0l,r(s, i^) has been brought under control, and the exploration 
of the resulting W{uj) has commenced. Further results are expected shortly. 



III. INTERACTING PROBLEM 

The transformation leading to (0) also applies to the general if in (P), and scaled param- 
eters (e. A, A) are readily written down in analogy to a, u. However, after this, no further 
exact simplification is apparent for H, so that the interacting problem is considerably richer 
than the uncoupled one. However, there is an alternative form 

H = Hs- Ina, + {Xz(T,+X^a^)s , (15) 

in which the free qubit Hamiltonian, not the interaction, is written in diagonal form. Specif- 
ically, defining cotr] = e/A, < i] < n, one has H = e*''°"f /^ife"""^"""/^, provided that 
Xz = XcosT] and A^,. = — Asinr^. Thus, the spectra of H and H are identical, and their 
eigenfunctions trivially related. 



A. Analytics 

The above, while seemingly trivial, has an intriguing consequence for the analytical study 
of 0, which now becomes a well-defined problem in resonance perturbation theory. Rou- 
tinely following the standard prescription "if = ifqb + ifreadout + Hmt \ one is led to expand 
in ifint, i.e., in A. However, Q is also "solvable" (in terms of the resonances of H^) for 
finite A and A = 0, in which case the spinor components decouple. By ()15p. we see that a 
A-expansion assumes A^^ — at arbitrary A^, while the A-expansion assumes both X^^X^ — > 0. 
Conceptually, H is appealing: it shows that the QND component A^ of the coupling can be 
readily handled, while the non-QND (henceforth QD) component A^^ is nontrivial. 

In contrast, the A-expansion, which mixes both components, places an unnecessary re- 
striction on the QND component, and obscures the fact that a QD perturbation (either A^ 
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or A) affects the resonance frequencies only in second order. ^ Still, we give for reference 



duj{e, A, A) 



COST] 



A=0 



(16) 



where we used du{Q,, Xz=0, X^) /dXx\x^=o = 0. For the qubit in its ground state, one has 
cos?7 = {az) and duj^fl, Xz, Xx=0)/dXz\x^=o = i^S' "with cjg as in (fT^. where quantities 
pertaining to the uncoupled SQUID will be henceforth denoted with an "S" subscript. 
For an excited qubit, both these relations have the opposite sign. Hence, in either case, 
one obtains Su = Xuj'^{az) for the first-order frequency shift (where on the rhs only a;g 
has an imaginary part): self-averaging {Hamiltonian- dominated ^) measurement is a direct 
consequence of the weak-coupling assumption. The SQUID's working point a strongly affects 
the measurement sensitivity through tUg. 

Conversely, let us consider the first-order effect of the qubit splitting Q in the presence 
of a finite coupling A. For only QND coupling, the qubit spin always has a definite value, 
which we denote by Sz for H and Sz for H, so that Sz = Sz sgn e = Sz sgn A^. One now has 



cu(e,A,A)b^o~cu(0,0,A) + 
= LJsia+Xsz) ■ 



9cu(0,0,A) ^ ^duj{0,0,X) 



de 



dA 



esz 
2 



(17) 



where on the first line, the second term is trivial since it is evaluated in the absence of 
tunneling, and the third term has been argued above to cancel. One can interpret ()17p 
by saying that strong coupling to the SQUID has collapsed the qubit in the pointer basis 
(az eigenbasis) and suppressed tunneling — detector- dominated measurement |5l]. 

For a systematic approach, one should relate and work out the A- and A2,.-expansions 
analogously to (fT^: 

d^uiQ cos 7], Q sin r/. A) 1^^^^^ = d'^u{Q, X cos r], -A sin rj) \^^^^^ 



9A2 



A=0 



XzSzUj'sia+XzSz) + X] 



dXl 



(18) 



A^=0 



That is, in general the two schemes yield equivalent information, but the A-expansion be- 
comes trivial for A = (no coupling in any order), while the A^^-expansion becomes trivial 
for = (QND in any order), with of course d^uj = in either case since the problem is 
77-independent unless both A and Q are nonzero. Note also that, from its prefactor in ()18p. 
one cannot conclude that d\uj\A=o diverges for e — 0, because in this limit the two terms 
on the last line of (fTHj) cancel each other. Rather, d\u!\A=e=o should be calculated directly 
from (PUI) below, and cannot be trivially related to other expansion coefficients. Analogous 
remarks apply to dl^uj\x^=\^=o. 

Next, the actual expansion should be carried out. This is sketched for the A-expansion 



of the resonance with = 1; other cases are analogous. In O(A0), V'} 



(0), 



A) = V's(«+A) 



^ Note that the effect of A^; or A must be even, since their sign is a matter of convention, or draw a picture 
on the Bloch sphere. 
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and ?/'|°'' = 0, while a;^°^(a,e, A) = ti;s(Q^+'^) ~ In ^{^)^ the couphng of the spinor 
components means that is a source term in the equation for which is solved as 

'ilj^^\s]a,e,X) = —^ JduG(s,u;uJs{a+X)—e]a—X)ilJs{u,a+X), (19) 

where the Green function (I12p of the free SQUID is simply denoted as G since there is no 

risk of confusion. Finally, ipf^'^ its turn figures as a source term in the equation for 
from which the frequency shift can be determined in much the same way as in (jl4p . The 
final results are 
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dXl 



_ f^dsdu'4js{s,a+Xs^) G{s,u]LJs{a+Xsz)-eSz]a-Xsz) iJs{u,a+Xsz) ^^^^ 
A=o 2 J^dsijl{s,a+Xsz) 

_ 2 J^ds du sepsis, a+XJz) G{s, u; uJs{a+XzSz)-^h] a-X^s^) wpsiu, a+XzSz) 
A,=o Jcdsi:l{s,a+XzSz) 

(21) 



There seems to be an interesting resonant enhancement ii (jj^\a+Xsz)—eSz ~ uj^\a—Xsz), 
where the SQUID initially is in mode j, while k denotes another mode. Namely, in this 
case G in ()2Up is evaluated near a pole. However, caution is needed, for this "resonance" 
condition could be met even if A = 0, where one knows that ()20|) in fact is trivial. 

Note that, if the free SQUID's eigenf unctions were complete, the Green function (or 
resolvent) would have a representation of the form G{s,u]U!) ~ '^^ips,n{s)ips,n{u)/ {uj—ujn)- 
In such a case, ()2()j) and ()2H1 would reduce to a more familiar ratio of squared matrix elements 
and "energy denominators" , summed over the mode index . 

It is instructive to verify the predicted equivalence (jl8|) directly from these integral ex- 
pressions. For that purpose, we write (0) as 

-ips{s, a)" + [ais - - LJi]tps{s, a) = [ujs{a) - ui + {ai-a)s]ips{s, a) , (22) 

for arbitrary ai,uji G M. One can invert the differential operator on the Ihs, yielding 

ijjs{s,a)= duG{s,u;uJi; ai)[iUi — LUsict) + {a—ai)u]^jjsiu, a) . (23) 



Jc 

Applying twice to (PT|). first to the w-integral and then to the s-integral, one verifies that 
PU)) and (jnH) are indeed related by (fTH|) (in particular, also when one of the two expressions 
becomes trivial). □ 

The above represents substantial formal progress, but further analysis is needed to eluci- 
date the physics contained in and (PT|) . Presumably, one should try to evaluate these, 
in particular their imaginary parts, in the semiclassical limit n <^ Ns, and ideally make a 
comparison with nonperturbative numerical results for the interacting system. 



B. Degenerate perturbation theory of the hnear-detector regime 

The above allowed for an arbitrary angle t] between the energy and pointer bases only 
in two disjoint cases: infinitesimal A as in (fTB|). or infinitesimal Q as in ()17p. One is led to 
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suspect that a unified treatment is possible if both X and Q are small, but this regime cannot 
be directly extracted from the above. Namely, in the leading order A = = 0, the problem 
reduces to the free SQUID, and hence is degenerate in the spin variable. One also notices 
that the central formulae (|2Up and (|2H) become singular in this limit, for the Green functions 
G are needed infinitesimally close to their resonance frequency, where they diverge. 

Therefore, here we give a separate calculation, expanding in e. A, and A simultaneously. 
Physically, this (presumably) means that both the qubit splitting and the typical interaction 
energy are very small compared to the SQUID plasma frequency. In its turn, this means that 
only one of its (noninteracting) resonances is needed to describe the SQUID (typically the 
longest-lived "ground-state" one), i.e., that the SQUID follows the qubit adiabatically — all 
in all, adiabatic elimination is just a form of perturbation expansion. In the first order, 
which will be considered here, this moreover agrees with the linear- detector regime, when 
the qubit perturbs the detector only weakly, without necessarily assuming the reverse. 

Typically, degenerate perturbation theory starts with considering the matrix elements of 
the perturbing operator in the degenerate subspace. Here, some care is needed since there 
is no concept of a basis; however, one can proceed in analogy with (jl4|) . and the product 
(•, •) will emerge naturally. Thus, consider the leading correction to (0) analogous to (fT^. 

[Hs - cosia)]\E) + [-5iu - i(ea,+Aa,) + Aa,s]|^) = , (24) 

where 6lj is the leading correction to the resonance frequency ujs{a), |S) = (^f^i)''^ 
corresponding correction to the lowest-order wavefunction |\E') = (^q^ips, and where p,q are 
as yet undetermined. As before, Hs — t^s(«) is not invertible, so a solution |S) satisfying 
both boundary conditions only exists in a special case, namely if 

-i£+A^ _-lA \ fp\_,Jp 



This is a key result: the resonances 6u! follow from a matrix equation, involving 
both the qubit parameters e, A and the complex linearized detector response cjg, as 
5uj = ±^^yJe^^^2XuJ^)^~+~^. Depending on the relative magnitudes of and A|u;g|, the 
Hamiltonian- and detector-dominated regimes are now trivially extracted, as will be sketched 
below for the purpose of interpretation. 

First, however, let us pursue the determination of Assuming that p, q, and 6uj are 
related as in (j25j) . one can rewrite (j24j) as 

[Hs - u;s{a)]\E) + A(s - {s))a,\^) = , (26) 

and comparison with (fT!^ yields the solution |S) = A(t^9q,|\E'). That is, to this order, the 
corrected wavefunction can be written as 

I*.) « I*) . IE) . (*(-^)) . (27) 

In agreement with the above adiabatic picture, (j27j) gives the SQUID wavefunction as the 
instantaneous ez^'enfunction corresponding to the qubit flux state. In particular, and un- 
surprisingly, the system's wavefunction is essentially always entangled, with the degree of 
entanglement characterized by the coupling A. The exception are eigenstates with p — >■ or 
g — 0, that is, QND or detector-dominated measurement of a classical flux state. 
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Let us discuss the two limiting cases. In the detector-dominated regime Q ^ ■^I'^sL 
one finds Suq ^ Xu'^ ~ 1^ + ^"^/^Xuj's and {p,q)Q oc (1,— A/4Aco'g — Ae/8(Acjg)^), while 
Sui = —Sujq and {p,q)i oc (A/4Aci;g + Ae/8(Aco'g)^, 1). That is, the qubit has "collapsed" in 
the pointer basis; the frequencies are mainly determined by the detector response and the 
bias e, while tunneling ~ A^ is strongly suppressed with a factor (Acjg)^^. Conversely, in the 
Hamiltonian-dominated regime Q ^ -^li^sl' ~ ^l^^ ~ Xu'^e/Q + [Xugj^A'^/fl^]. 

The first two terms are ±|f2 + Aa;g(cr^), familiar from below ()16|) . Rewriting the last term 
as 5a; = • ■ ■ ± {XxUj'gY/fl, we tentatively interpret as follows: for Wg mostly real, A^ causes 
SQUID-assisted tunneling, increasing the effective tunnel splitting. For Wg mostly imaginary, 
however, the measurement of may be said to cause a Zeno effect, reducing the effective 
tunnel splitting. The corresponding wavefunctions are readily written down when needed. 

As clearly as (j25p shows the transition from iJ-dominated to detector-dominated mea- 
surement, it does not address all conceptual questions. Namely, it always describes weak 
measurement; on the opposite side, one has single-shot readout. "Single-shot" means that 
for a suitable waiting time, decay for one state is almost complete, while for the other state 
it is negligible, i.e., |Ima;o| In general, one presumably should introduce a dis- 

criminating power based on the ratio of these decay rates. In this framework, one can study 
issues such as whether single-shot measurement can be if-dominated, and the existence of 
a dimensionless quantum efficiency < ?7 < 1 like for continuous measurement. 



C. Numerics 

Since the uncoupled Schrodinger equation for Hs already is implemented numerically 
as two coupled first-order equations, the integration of should only involve a trivial 
generalization. Also, the cubic potential term continues to dominate asymptotically in the 
presence of a qubit, so the boundary conditions and choice of integration contour C carry 
over unchanged. The matching, however, needs more care. Based on the trivial cases 
A = or A = 0, one now expects two linearly independent solutions for each of and 
I^'r) satisfying the corresponding boundary conditions. During the numerical integration, 
we have no a priori knowledge which solution is the correct one, so we have to integrate 
both, and determine a fundamental system |\1'l'*)} of independent spinors, obeying 

I^L H-^)) — * for s — > — oo (j = 1,2), and analogously for L ^ R. Subsequently, the question 
becomes for which uj is there a superposition of the I^I'l '') matching onto a superposition of 
the of course for both spinor components and for the first derivatives as well as the 

function values. In other words, one has to find the zeroes of the Wronski determinant 

V'}l(«o,^) ^^iso,^) -iPriiso^uj) iP\'^tI{so,uj) 

^{l(so,^) V'{?(so,^) ^[r(so,^) ^{JJ(so,a;) 

iIj\^^'{so,uj) ip\f^'{so,uj) iIj\]1'{so,uj) V'}JJ'(so,a^) 

V'{l'(so,c^) ^Jl'(*o,c^) ip^^'{so,uj) iIj^^'{so,uj) 

Ignoring the spin index, the form (|7j) for the uncoupled Wronskian is seen to be precisely 
the ( (3'i) (3 3)) minor determinant of (j2Hl)- 



(28) 
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APPENDIX A: TOY MODEL 

Consider a particle with mass m = | on the positive a;-axis, 

+ U{x)ij{x) = ujip{x) . (Al) 

The origin a; = is supposed to be a hard wall, yielding the boundary condition tp{0) = 0. 
For X > 0, the potential is a finite well, U{x) = V6{a—x), with 6 the unit step function and 
V < < a. Depending on V, a, the problem has a finite number (possibly zero) of bound 
states with < a; < 0. These are to be determined by the Wronskian method of ((Tj), with 
x = a being the obvious choice for the matching point. On < a; < a, the "left" function 
ipi^(x,uj) = sm{\/uj—Vx)/\/uj—V is readily determined from the nodal condition. The same 
is true for the "right" function = e~^'^~^^ on x > a, but only because the closed- 

form solution enables us to deal with the singularity at a; = +oo; instead, one can examine 
how "inward" integration works in this example. First, pick some trial initial values, say 
ipnixaju) = —iP'y^{xq,u) = 1; in general, {ipnixcu), iIj'^^{xq,lu)} could be left undetermined. 
Second, solve the Schrodinger equation ()A1|) on x > a with these trial initial values, e.g. 
normalizing the result such that tpji{a,uj) = 1. [In fact, as usual, the logarithmic derivative 
il)'^^{a,Lj)/ipii{a,uj) is all that matters for Wronskian matching.] For the toy model, this can 
be done analytically, while in general it would be a numerical task. Third, consider ip^i as 
xq becomes large, while the trial initial values {iPr{xo,uj), ip'j^{xo,uj)} are kept fixed. One 
observes that the error committed by choosing a finite xq decreases exponentially. 

In slightly more technical language, the above shows that, even though x = cxd is an 
irregular singular point of the Schrodinger equation for a free particle, one can still effectively 
implement the boundary condition '?/'Pi,(x=+oo, Ci;<0) = numerically. 

To complete this exercise, one could now perform the matching and calculate approximate 
eigenvalues {(Dj(xo)}, studying their convergence to the true spectrum {uj} as xq — > oo. 
Here, ()A1|) is instead pursued analytically, as a toy model for resonance states of the 
Schrodinger equation. For < as above, this is illuminating, since the bound states 
now feature in both types of spectrum of the system: they are the discrete part of the con- 
ventional real hermitian-operator spectrum, and the undamped members of the resonance 
spectrum. For (Re) a; > 0, the continuous conventional spectrum and the damped discrete 
resonances offer complementary descriptions, and correspond to different definitions of an 
eigenvalue. Moreover, the model can be considered for V > as well, in which case there are 
also resonances. This clarifies that the resonances are caused by wave scattering rather than 
by classical confinement: in quantum mechanics, confinement is not absolute at energies 
below the pertinent barrier energy (due to tunneling) nor is it absent at higher energies, for 
which classical particles escape unimpeded. 

Performing a scale transformation as in the main text, we can take a = tt w.l.o.g., upon 
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TABLE I: Numerically determined roots of the Wronskian in ()A2|1 . corresponding to the first few 
resonance frequencies of the square- well model ()A1() . 



which the Wronskian becomes 



W{uj) 



sm 



cos 



(Vcu-Vtt) 



(A2) 



In order for this expression to properly account for the decaying boundary condition for 
X — s> 00 and C(j < 0, one has to define ^/uJ with the branch cut on the negative imaginary 



e.g., Mathematica, this can be done as 



Sqrt[ 



-luo 
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The real roots 



of W{ijj) can be found graphically as in elementary quantum mechanics; one finds that 
there are n bound states for — i(2n+l)^ < V < — |(2n— 1)^. In general, it is convenient 
to first explore the spectrum by plotting |Vr(a;)|~^; once the resonances have been located 
approximately, numerical minimization of |Vr|^ reliably determines their accurate value. The 
first few resonances for two instances of the model are given in Table HI Based on further 
exploration, it is conjectured that Reci;„ > in all cases. 

Clearly, the low-lying bound states for large negative V can be found by approximating 



the well as a box with hard walls, as Un 



n 



+ V . Here, we focus on the less-known 



high-order modes, with |a;| ^ \V\. For this purpose, set u = uj — V and expand in V. The 
resonances follow from e^*'^^ ~ 4z//V^ + 1, which is asymptotically solved by 



2 -"^ 1 
z/^ ~ m — I — m 



Am/ 



\t? m 



+ 



(A3) 



For < 0, different choices of the logarithm simply yield predictions for different modes. 
Comparison with Table H] shows good agreement for either choice of sgn V, although the 
presence of various logarithmic terms in ()A3|) means that convergence is slow. Interestingly, 
the dominant term is the same as for an infinite well. For V ^ 0, the decay rates 
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Since W{uj) in (jA2|) is even in \/uj—V, there is no branch ambiguity for the latter. 
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diverge as expected, but perhaps counterintuitively, for any fixed V the quahty factor Qm = 
Reci;m/2|Imci;m| increases without bound. 

For comparison, we briefly mention the model U{x) = K6{x — tt), for which the modes 
are asymptotically 

2 .n . f 2in\ n , . ^, 

Comparing with ()A3|1 . one sees that the stronger potential singularity (J-peak versus step) 
causes stronger scattering, and hence a decay rate which asymptotically is a factor two 
smaller. 

A closed-form Wronskian (in terms of hypergeometric and related functions) should 
also be obtainable for many models with a potential consisting of two constant, linear, 
or (inverted-) harmonic pieces, on either the full or the positive x-axis. However, unlike any 
of these, the cubic potential is infinitely smooth, so that its asymptotic scattering properties 
(reflected in the damping rates of high-order resonances) are expected to be qualitatively 
different. 
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